Search for long lived heaviest nuclei beyond the valley of stability 
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The existence of long lived superheavy nuclei (SHN) is controlled mainly by spontaneous fission 
and a-decay processes. According to microscopic nuclear theory, spherical shell effects at Z=114, 
120, 126 and N=184 provide the extra stability to such SHN to have long enough lifetime to be 
observed. To investigate whether the so-called "stability island" could really exist around the above 
Z, N values, the a-decay half lives along with the spontaneous fission and /3-decay half lives of such 
nuclei are studied. The a-decay half lives of SHN with Z=102-120 are calculated in a quantum 
tunneling model with DDM3Y effective nuclear interaction using Q a values from three different 
mass formulae prescribed by Koura, Uno, Tachibana, Yamada (KUTY), Myers, Swiatecki (MS) and 
Muntian, Hofmann, Patyk, Sobiczewski (MMM). Calculation of spontaneous fission (SF) half lives 
for the same SHN are carried out using a phenomenological formula and compared with SF half 
lives predicted by Smolanczuk et al. Possible source of discrepancy between the calculated a-decay 
half lives of some nuclei and the experimental data of GSI, JINR-FLNR, RIKEN are discussed. In 
the region of Z=106-108 with N~ 160-164, the /3-stable SHN fo6-S'ffi62 is predicted to have highest 
a-decay half life (T a ~ 3.2/irs) using Q a value from MMM. Interestingly, it is much greater than the 
recently measured T a (~ 22s) of deformed doubly magic io|^ s i62 nucleus. A few fission-survived 
long-lived SHN which are either /3-stable or having large /3-decay half lives are predicted to exist 
near 294 110i84, 293 110is3, 296 112i 8 4 and 298 114i 8 4. These nuclei might decay predominantly through 
a-particle emission. 
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PACS numbers: 27.90.+b, 23.60.+e, 21.10.Hw, 21.30.Fe 



I. INTRODUCTION 

The theoretical studies of properties of heaviest nuclei 
during the past few decades have drawn considerable 
attention of experimentalists to investigate the existence 
of superheavy nuclei beyond the valley of stability. 
Since the macroscopic [H,IH description on the basis of 
liquid drop model (LDM) does not take the shell effect 
into account, it fails to explain the variation of fission 
barrier height of heavy nuclei with the increase of the 
fissility parameter (~ Z 2 /A). However, according to 
modern nuclear theory, hindrance to the fissioning of 
heavy nuclei would be enhanced due to the presence of 
deformed and spherical shell closures. Different semi- 
microscopic approaches e.g. macroscopic-microscopic 
model (MMM) @, H i, H and its modification Q 
include pairing and nuclear shell effects 0] to reproduce 
the properties of ground and deformed states of nuclei. 
Many purely microscopic [H, [n| EL G2, EH descriptions 
like Hartree-Fock-Bogoliubov (HFB) model with zero 
range forces of Skyrme type [lil . [l5j or finite range 
forces of Gogny type [H, E3] and relativistic mean field 
(RMF) [la, [l9( theory predict the possible deformed and 
spherical neutron shell closures at N=162 and N=184 
2C | respectively. Since strong influence of nuclear shells 
2l| in the region of superheavy elements might make 
sufficiently long lived SHN to be observed, the search 
for heavier elements in the natural samples was started 
[H, H [H, HI about thirty years ago. 

Experimental investigations in finding the SHN around 



Z=107-118 have been pursued mainly at three different 
places: Gesellschaft fur Schwerionenforschung (GSI) 
in Darmstadt (Germany), Joint Institute for Nuclear 
research (JINR) in Dubna (Russia), and RIKEN, Japan. 
In the beginning of the 1980's the first observations of 
the elements with Z= 107-109 were made at GSI [26j. In 
1994, a-decay chains were observed from nucleus 269 110 
[27] | and later on a-d ecay chains from nuclides 271 110, 
272 111, 277 112 mEiUjl, 283112 S3 were detected at 
GSI. 



While RIKEN claimed discovery of the 278 113 
SHN [H, HH, it also reconfirmed the a decay 
chains from 271 110 0, 272 111 H and 277 112 
[36| . Observations of the a decay chains of nuclei 

294 -Qg 290-293-Qg 288,287 -Qg 286-289 -q 4 
285,283]_]_2 [32J ' 278-280 -qi 
274-276 10 Q 275 JOS 272,270 ^ 

by jinr go, EH, El] • 



282-284 113) 

273 ' 28i no n m, 

271 106 were reported 



Recently, the 270 Hs (Z= 108, N = 162) SHN has been 
produced in the 26 Mg + 248 Cm reaction j44|. According 
to the theoretical calculations 0, Hf|, nucleus 27t) Hsi e2 
(Z=108) should have the features of "deformed doubly 
magic" nucleus. Most of the heaviest nuclei are expected 
to be deformed due to partial filling of large nuclear 
shells by outer nucleons. Dvorak et al. measured the 
energy (E a ) of a particle emitted from 270 i/si62 and 
used the value of E a to calculate Q a (9.02 ± 0.03 MeV) 
for the a decay of 270 i?si62- A phenomenological 
formula [46|] estimated the a decay half- life (~ 22 s). 
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Earlier, it was believed [H, H3, IH, El| that traditional 
spherical superheavy nuclei might form an "island of 
stability" centered around 298 114ig4 separated from 
the "peninsula" of known nuclei by a region of deep 
instability. Due to both deformed neutron shell and 
proton shell effects at Z=108 and N=162 the extension 
of the peninsula of known nuclei might connect the 
stability island of spherical superheavy nuclei around 
doubly magic spherical Z=114 proton shell and spherical 
N=184 neutron shell. Since fission barrier and shell 
effect play very important role for the existence of long 
lived superheavy nuclei it is crucial to determine the 
fission barrier and half life of fissioning nucleus with a 
good accuracy. It is well known that very small barrier 
height against fission can break the nucleus into two 
fragments immediately after it is formed. The a. decay of 
superheavy nuclei [H(| HH [H, [H, [Hi], [HH], H(| is possible 
if the shell effect supplies the extra binding energy and 
increases the barrier height of fission. /3-stable nuclei 
having relatively longer half life for spontaneous fission 
than that for a decay indicates that dominant decay 
mode for such SHN might be a-decay. 

In our previous works [HI], [H, [H3, [H, we showed 
the applicability of the microscopic calculation in 
predicting the a decay half lives of SHN from a direct 
comparison with the experimental data [4(| Hi], [42l |. 
However, as a number of SHN were predicted to have 
relatively large a decay half lives, it is necessary to 
find out whether those SHN would survive the fission 
[Hf^ and /3-decay. In such cases those nuclei can be 
detected in the laboratory through a decay. This work 
explores the possibility of finding long lived SHN by 
comparing the calculated a decay half lives (T a ) with 
available theoretical spontaneous fission (SF) half lives 
[HH [HU, calculated /3-decay half- lives (Tp) (63[ and the 
experimental data on SF. The a decay half lives of SHN 
with Z= 102-120 are calculated in a quantum tunneling 
model with DDM3Y effective nuclear interaction using 
Q a values from three different mass formulae. 

A brief outline of the methodology of the present cal- 
culation is presented in section II. S pon taneous fission 
half lives from both phenomenological [64l IHH| and micro- 
scopic approach [HI], [H2] are given in section III. Finally, 
in section IV, results and discussions and in section V, 
summary and conclusion are presented. 

II. FORMALISM 

The a decay half lives are calculated in the frame work 
of quantum mechanical tunneling of an a particle from 
a parent nucleus [H3|- The details of calculation of the 
a decay half lives of superheavy nuclei were described 
in our earlier works [H?], [H, [H§|. The required nuclear 
interaction potentials are calculated by double folding 
the density distribution functions of the a particle and 
the daughter nucleus with density dependent M3Y ef- 
fective interaction. The microscopic a-nucleus potential 



thus obtained, along with the Coulomb interaction po- 
tential and the minimum centrifugal barrier required for 
the spin-parity conservation, form the potential barrier. 
The spin-parity conservation condition in a decay process 
is fulfilled if and only if 

J=J X +J 2 + 1, 7r = 7ri.7ra.(-l)', (1) 

where J, Ji and J 2 are the spins of the parent, daugh- 
ter and emitted nuclei respectively, tt, m and 1^2 are the 
parities of the parent, daughter and emitted nuclei re- 
spectively, and 1 is the orbital angular momentum carried 
away in the process. This conservation law, thus, forces 
a minimum angular momentum to be carried away in the 
decay process. Consequently, contribution of the angular 
momentum gives rise to a centrifugal barrier 

Vi = h 2 l(l + \)/{2 i iR 2 ) (2) 

where /i is the reduced mass of the daughter and emitted 
nuclei system and R is the distance between them. 

The half lives of a disintegration processes are 
calculated using the WKB approximation for barrier 
penetrability. Spherical charge distributions have been 
used for calculating the Coulomb interaction potentials. 
Most of the experimental Q-values of a decay (Q a ) are 
obtained from experiments done in GSI, Germany and 
JINR, Dubna. For theoretical Q-values, mass formulae 
from KUTY [H(|, Myers-Swiatecki [H3| and Muntian et 
al. @, [HI, 11] are used. 

The experimental decay Q values (Q ex ) have been ob- 
tained from the measured a particle kinetic energies E a 
using the following expression 

Qe X = {- r ^-)E a + (65.3Z p 7 / 5 - 8O.OZ2/ 5 ) x lO" 6 MeV 
A p — 4 

(3) 

where the first term is the standard recoil correction and 
the second term is an electron shielding correction in a 
systematic manner as suggested by Perlman and Ras- 
mussen [70| . A p and Z p are mass number and atomic 
number of parent nucleus respectively. 

The theoretical decay Q values Qth have been ob- 
tained from theoretical estimates for the atomic mass 
excesses 0, [HI, [H3: [H3, HH using the following relationship 

Q th = M - (M a + M d ) = AM - (AM a + AM d ) (4) 

which if positive allows the decay, where M, M a , Md 
and AM, AM a , AMd are the atomic masses and the 
atomic mass excesses of the parent nucleus, the emitted 
a particle and the residual daughter nucleus, respectively, 
all expressed in the units of energy. As Q^-value appears 
inside the exponential integral as well as in denominator 
of the expression [Hj, [58[ of a decay half lives, the entire 
calculation is very sensitive to the Q-values. 
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III. PHENOMENOLOGICAL AND 
MICROSCOPIC CALCULATIONS FOR 
SPONTANEOUS FISSION HALF LIVES 



Spontaneous fission of heavy nuclei was first observed 
by Flerov and Petrjak in f940 [7l| from 238 U nucleus. 
Spontaneous fission and a decay are the main decay 
modes [lH, [z!, [ill °f superheavy nuclei. /3-decay could 
be another possible decay mode for the superheavies 
lying beyond the /^-stability line of the nuclear chart. 
However, since the /3-decay proceeds via weak interac- 
tion, the process is slow and less favored and energy 
involved (released) is also less compared to spontaneous 
fission and a decay which proceed quickly via strong 
interaction making these processes more probable. For 
heaviest nuclei, the mutual repulsion of electric charge 
is higher than surface energy of the nucleus arising 
from the short range nuclear forces. On the basis of 
liquid drop model, Bohr and Wheeler [75| described the 
mechanism of nuclear fission and established a limit 
for Z 2 1 A ~ 48 for spontaneous fission. Beyond this 
limit nuclei are unstable against spontaneous fission. 
According to this fact, no nucleus beyond Z ~ fOO can 
exist due to very small fission barrier. Both theoretical 
and experimental investigations on superheavy nuclei 
(SHN) support the fact that a bound superheavy nucleus 
can be formed only due to shell effects. 



ing formula (65| : 



logw(T 1/2 /yr) = 21.08 + a 
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where c x = -548.825021, c 2 = -5.359139, 
c 3 = 0.767379, c 4 = -4.282220, v = for the 
spontaneous fission of even-even nuclei and v = 2 for 
odd A and odd-odd nuclei. The seniority number v was 
introduced for taking the blocking effect of unpaired 
nucleon on the transfer of many nucleon-pairs during the 
fission process. A variation of SF half lives calculated 
by this formula (eqn.5) with increasing neutron number 
for different elements from Z=102 to Z=120 are shown 
in Fig. 4 and Fig. 5 for comparison with the a decay 
and SF half lives predicted by our calculation and by 
Smolanczuk et al. [6ll |62| respectively. Moreover, in 
this work, spontaneous fission half-lives calculated in 
a dynamical approach using macroscopic-microscopic 
method (MMM) of Ref.[6l|, Isl are also used (Figs. 2-5) 
to find out the "island" where the predicted SHE could 
survive the fission. 



A simple semi-empirical formula on spontaneous fission 
half-lives for even-even, odd A and odd-odd nuclei in the 
gro und state was proposed by W.J. Swiatecki in 1955 
76]. By including the deviation of experimental ground 
state masses from a smooth reference surface based on 
liquid drop model, Swiatecki successfully reproduced the 
experimental data with his semi-empirical formula [76| . 
Microscopic calculation of spontaneous fission half-lives 
is very difficult due to both the complexity of the fission 
process and the uncertainty of the height and shape of 
the fission barrier [77j . Ren and Xu (64|, [65[ generalized 
the formulae of spontaneous fission half-lives of even-even 
nuclei in their ground state to both, the case of odd nu- 
clei and the case of fission isomers [65| . The spontaneous 
fission half-lives of odd-A nuclei and of odd-odd nuclei 
in the ground state were calculated by using the gen- 
eralized form of the Swiatecki's formula and by a new 
formula where the blocking effect of unpaired nucleon on 
the half-lives were taken into account. By introducing a 
blocking factor or a generalized seniority in the formu- 
lae of the half- lives of even-even nuclei, the experimental 
fission half-lives of odd-A nuclei and of odd-odd nuclei 
with Z=90 to Z=108 were reasonably reproduced with 
the same parameters used in ground state of even-even 
nuclei. 

In the present work, spontaneous fission half-lives for 
both neutron deficient and neutron rich isotopes of ele- 
ments Z=102-120 have been calculated using the follow- 



IV. RESULTS AND DISCUSSION 

In this paper, our main aim is to check whether the 
predicted SHN do really survive against the spontaneous 
fission (SF). Gamow- Teller /3-decay half lives (Tp in 
figs. 2-3) obtained from a microscopic quasi-particle 
random phase approximation with single particle levels 
by Moller, Nix, Kratz [63| has also been used in this 
work to check the possibility of /3-decay of SHN. In 
earlier works [5l|, [3 [S3, [H, H3|, we were able to well 
reproduce the available experimentally measured a 
decay half lives following the semi-classical quantum 
tunneling method with double folded density dependent 
effective M3Y interaction. In Ref. [58j, T a of about 
314 nuclei werepredicted using Q a values from MMM 
(Qm) 0> HI, [6{| and modified liquid droplet model 
of Myers-Swiatecki [67| to show the necessity of more 
accurate mass formula in determining the Qc-values 
with a good accuracy at least correct upto lOKeV for 
heaviest nuclei. But the spontaneous fission survivability 
of those nuclei was not checked in our previous work. 
For this reason, although some SHN (Z=102-120) are 
long-lived against a-decay [StI [Hj], they may not live 
long if they have shorter SF half lives. For example, at 
N=162 the magnitudes of T a using Qm for the elements 
No (Z=102), Rf (Z=104), Sg (Z=106), Hs (Z=108) are 
3.06 x 10 9 s, 4.10 x 10 6 s, 1.16 x 10 4 s, 28s respectively 
[58| . If SF half lives (Tsf) of such nuclei are about the 
same order or relatively longer than their respective T a , 
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only then a significant fraction of nuclei can survive 
fission and decay by a emission. In such cases a-decay 
chain followed by SF of one of the product nuclei may 
be observed. Therefore, accurate estimation of Tsf of 
SHN is essential to know whether predicted long-lived 
SHN against a emission do really exist against SF. 

Gupta and Burrows [78| summarised the measured 
ground state values of spontaneous fission and a decay 
half lives with Q a for heaviest nuclei having mass 
number A=266-294. These values were taken from ex- 
perimental measurements carried out at GSI, Germany, 
and JINR, Dubna. In Table I, a comparison between 
experimentally measured and calculated half lives using 
DDM3Y effective nucleon-nucleon (NN) interaction in 
a WKB framework is shown for sixteen superheavy 
nuclei. Most of these nuclei (twelve out of sixteen) were 
not addressed in our earlier works [5l|, [H, [57], [H, l59j . 
In Fig.l the theoretical Q Q -values from different mass 
formulae prescribed by MS [67j], KUTY[66j], Muntian 
et al. (MMM) 0, [H, [1| are employed to calculate 
a-decay half-lives of the same sixteen SHN presented in 
Table I. The measured Qa-values are used to calculate 
a decay half lives {T^ M3Y [Q ex }) in Table I. Half lives 
calculated in the present work agree reasonably well with 
experimental data for most of the superheavies. Nuclei 
for which spin-parities of parent and daughter nuclei are 
not known, zero angular momenta (£ — 0) transfers are 
used. Since £ — gives minimum centrifugal barrier, 
probability of a-tunneling increases and consequently 
half life decreases. But, for some nuclei like 277 112, 
273 110, 269 110, 266 106 etc. angular momenta carried by 
a particles may not be zero. For these nuclei, calculated 
half lives are too small compared to measured values. 
In particular, for 277 112 calculated value (~ 0.093 ms) 
is less than the measured one (0.69 ms) by one order 
of magnitude. This discrepancy can be removed only 
if one can assume non-zero orbital angular momentum 
transfer in the a decay process. In our calculation, 
£ = 4 gives the better agreement of predicted half life 
for 277 112 nucleus with measured value. It is therefore 
important to determine the proper spin-parity of parent 
and daughter nuclei to enable correct centrifugal barrier 
for a decay half life calculation. Although, higher 
i'-values are shown for some nuclei for better agreement 
with the experimental result, but more experiments on 
the same nuclei (shown in table I) with higher statistics 
are needed for reconfirmation of the data and measured 
a decay energies since reaction cross-sections are of the 
order of pico-barns (pb). 

In this work, existing SF calculation with a micro- 
scopic approach [6ll |62| | have been used to find the region 
of long lived fission survived nuclei in the SHN region 
of the nuclear chart. In Fig. 2 and Fig. 3, comparisons 
between calculated T a , Tp and Tsf are shown only for 
even-Z elements with proton number Z=102-120 since 
T a and Tsf calculations of Ref. [U [12] are valid only for 



even-even nuclei. A few available observed data for both 
T a and Tsf are also shown for most of the elements in 
the plots of Figs 2-3. For Z=104, the highest value of 
T a (~ 4.1 x 10 6 s) according to calculation of this work 
using Qm appears around N=162 where much smaller 
value of Tsf calculated in Ref. [6l[ (~ 23 s) makes 
this nucleus 266 104i62 unstable against SF. Therefore, if 
synthesis of this nucleus is possible in the present-day 
setup, substantial fraction of 266 104i 6 2 would undergo 
SF within a few seconds. For Z=102, no such calculation 
on SF is available in Ref. [6ll. l62| . Since calculated Tsf 
and Qm values both are based on MMM, a decay half 
life calculations using DDM3Y effective interaction with 
Qm (T a M3Y Q-M) is preferably chosen to compare 
with Tsf ■ In addition to that, calculation within the 
same framework using Qkuty have also been presented 
in Figs 1-3. 

In case of Sg (Z=106) isotopes, DDM3Y with Q M 
predicts that the longest T a (~ 1.16 x 10 4 s ~ 3.2 hr) 
would be at N=162. It is comparable to Tsf (3.5 hr) 
of Ref. 61]. Therefore, a decay channel is one of the 
dominant decay mode of this 268 S , gi62 nucleus. If it is 
produced in the laboratory, it may be observed only for 
few hours (lifetime ~ 1.68 hr) since both SF and a 
decay half lives are small. It is to be noted that in the 
present work, lifetime of some SHN (either /3-stable or 
have large T ( a)are predicted by considering SF and a 
decay half lives only. Incidentally, 264 A0162, 266 Rj 162, 
26S Sgie2, and 270 Hsie,p are known to be either /3-stable 
or have very large Tp [63|]. Q a values using KUTY mass 
formula are not always very much reliable since it does 
not reproduce all the observed Q values with a good 
accuracy. But, this mass formula can be used to locate 
the region of possible existence of long-lived SHN where 
Q values from other mass formula are not available. 
It may be pointed out that the use of Qkuty hi our 
calculation shows reasonable agreement for several nuclei 
in Fig. 1. 

As Qm values are not available for more neutron 
rich isotopes of Sg, using Qkuty in present calcu- 
lation predicts T a ~ 3.71 x 10 15 s ~ 1.18 x 10 8 yrs 



at N=184. 



The SF half life of 

113 



2'«> 



Sg 



184 



is less 



(T SF ~ 4.07 x 10" s ~ 1.3 x 10 b yrs) than T a . Although 
Tp of this nucleus is not specified but it is expected to 
be large in ref.[(5j|- Hence, if synthesized, 290 Sgi$4 is ex- 
pected to have long enough lifetime (r ~ 1.27 x 10 6 yrs) 
to be observed in the laboratory. But it is still smaller 
than the age of earth ~ 4.5 x 10 9 yrs by three orders of 
magnitude. Calculated T a using Qm (28 s) for 270 -ffsi62 
is less than Tsf (~ 1-8 hr) and therefore a decay chain of 
such nucleus is expected to be observed. This prediction 
from our calculation is in good agreement with recently 
observed [44| doubly magic deformed nucleus 270 i?si 62 . 
Present calculation using experimental Q-value gives 
T a - 9.53s (Table I). 
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TABLE I: Comparisons between observed and theoretical (this work) Q-decay half lives using measured Q a . The experimental 
a-decay half lives {T^m P ) are taken from Ref.[78l] except the values with single (*) and double (**) asterisk symbols which are 
taken from Ref. [3l| and Ref. [44[ respectively. 



Parent 


Exptl Q- value 


Half-lives (Exp) 


This work(£ = 0) 


This work(^ 7^ 0) 




Ay 




rpEXP 

1 l/2 


rpDDMZY\ n l 
± l/2 LVea-J 


rpDDM3Y\s-\ 1 
J l/2 |>eea:J 


f 


283 112 


9.704 ± 0.015 


*6.9l 6 2 . 9 s 


4.67l 9 ; 4 4 9 s 






277 112 


11.594 ±0.055 


0.69^°'| 9 ms 


93^23 M s 


0.59±g ;^ ms 


4 


272 ln 


11.150 ±0.035 


3.8±l A s ms 


1.3ll°; 27 ms 






273 110 


11.37 ±0.05 


O.ntVol ms 


-7C + 21 
•0-17 M s 


0.13±r 3 ™ 


2 


271 110 


10.899 ± 0.020 


1.63l 9 ;^ ms 


930lgJ° flS 


1. 151^2 ms 


1 


270 110 


11.20 ± 0.05 


0.10i 9 ; l ms 


0.083l^ 2 g ms 


0.10l° ; 9 2 ^ 


1 


269 110 


11.58 ±0.07 


179l 2 f »s 


28lg 3 


88±H fis 


3 


267 110 


12.28 ± 0.11 


2.81S 3 ^ 








268 109 


10.486 ± 0.035 


2l±% ms 


i2.3±?,:2 ms 






266 109 


10.996 ± 0.025 


1.7±i;| ms 


750+H Ms 


1.3to.i ms 


2 


270 108 


9 qo+ 07 


3.6i?; 8 s 


i-36lS:S s 






**270 ^Qg 


9.02 ±0.03 


**22 s 


9-53+1 .86 s 






269 108 


9.315 ±0.022 


^■7-3.3 s 


3 1Q+ - 52 s 






267 108 


9.978 ± 0.020 


58t 23 ms 


45.5+g; 8 ms 






266 108 


10.336 ± 0.020 


2.3± ;| ms 


2.24l°; 27 ms 






267 107 


8.96 ±0.30 


17± 6 4 » 








266 106 


8.88 ±0.03 


&2+_\f s 


4.89l ; 2 ° s 


16t 4 3 s 


3 



From Figs 2-3 it is seen that the extra stability effect 
of neutron shell at N=162 almost disappears with the 
increase of atomic number and T a becomes of the order 
of millisecond to microsecond for the elements having 
Z > 110. On the other hand, in more neutron rich side 
around N=184 of elements Z=110, 112, 114 theoretical 
T a using Qkuty in present calculation are of the order 
of 10 10 s, 10 8 s, 10 6 s respectively which are much less 
than theoretical Tsf of the order of 10 12 s, 10 13 s, 10 13 s 
respectively. It must be noted that 296 112i 84 , 298 114i 84 
are /3-stable whereas 294 110i84 is predicted to have large 



Tp in ref. [63] due to its very small positive Q^-value. 
/3-stable nuclei and those with very large Tp are not 
shown in figs 2-3. 

Beyond Z=114 peak-value of T a plot (Fig. 3) around 
N=184 suddenly reduces showing a possible signature 
of spherical proton shell at Z=114. According to the 
present calculation, 298 114i 8 4 the so-called doubly 
magic spherical superheavy nucleus predicted by nuclear 
structure theory in the mid of 1960, has T a values of 
the order of 10 6 s and 5 x 10 2 s using Qkuty and Qm 
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respectively which are much less than Tsf ~ 10 13 s. 

In Figures 4-5, three curves describing spontaneous 
fission and a-decay half lives are shown in each graph 
excepting Z=102 for which calculation of spontaneous 
fission half lives in Ref. [6lL [62J , represented by Tsf SM in 
figs. 4-5, are not available. For the isotopes of elements 
Z=106, 108 the values of Tsf SM become of the order 
of one millisecond near N=170. However, the values 
of Tsf SM are comparable to the calculated a decay 
half lives of this work (T a M3Y Q-M) around N=154 
to 164 for Sg. In case of Hassium isotopes (Z=108, 
N=156-163), Tsf SM > T a M3Y Q-M reveals the fact 
that isotopes of Hs within this range may undergo a 
decay with a longest half life of the order of few seconds 
(~10-30 s). This is in good agreement with the recent 
experimental observation of a decay from the nucleus 
270 £/si62- But spontaneous fission half lives calculated 
using Eqn.(5) [ij, [65[ fall rapidly around N=160-170 for 
elements having Z=102-112. In fig. 5, SF half life values 
for neutron rich isotopes of elements having Z=114, 116, 
118, 120 are less than 10~ 15 s near N=180. For Z=114 
only, this calculation matches with microscopic results 
for some isotopes having N= 160- 170. This calculation 
contradicts the following two facts: (i) Experimentally 
measured value of T a for 283 112i7i is ~ 6.9 s (see 
Table I) whereas, SF half life for this nucleus calculated 
by using equation (5) is extremely low (— 10~ n s) 
indicating immediate fissioning of the nucleus, (ii) Since 
Tsf by Ren and Xu rapidly falls around N=160-170 
(fig.4) and N=180 (fig.5) for Z=102-112 and Z=114-120 
respectively, it could not explain the predicted extra 
shell stability of SHN at N=184 for heavier elements. 

From the present calculation of T a using Qkuty 
it seems that longer T a might be observed for more 
neutron rich side (Z > 116, N > 190) due to possible 
existence of neutron shell closure. On the contrary, from 
the trend of SF plots for Z=116, 118, 120 it appears 
that lowering of T$ f might destroy the existence of such 
neutron rich (N > 190) superheavy isotopes of elements 
Z=116, 118 and 120. Hence, the presence of long-lived 
SHN with neutron shell closure beyond N=184 may be 
ruled out. However, it is an important task to determine 
the SF and a decay half lives for N > 184 to confirm 
whether there is any possibility of neutron shell closure 
beyond N=184 for heaviest elements (Z > 116). 

In two graphs of Fig. 6, using Qm values in this cal- 
culation, variation of a decay with neutron number for 
both odd-Z and even-Z elements having Z=102-120 are 
shown. Plots of both graphs (a) and (b) of Fig. 6 clearly 
show peaks of T^-values around N=162 and N=184 for 
all elements with Z=102-120, which possibly indicates 
the neutron shell closure at N=162 and 184. This is in 
good agreement with the present-day knowledge of mi- 
croscopic theory. 



V. SUMMARY AND CONCLUSION 

The natural existence of superheavy nuclei is limited 
primarily by spontaneous a decay and spontaneous 
fission processes. A SHN in spite of having longer a 
decay half lives may undergo immediate spontaneous 
fission if the latter has a low half life. On the other 
hand, SF stable SHN may have shorter (~ 1/is or less) 
a decay half lives. In both the cases, such SHN may not 
be observed even if they are synthesized in the present 
day laboratory setup. The main aim of this work is to 
find out the fission-survived long lived SHN. In fact, if 
SHN have high degree of stability against both a decay 
and SF, we would be able to observe them if produced 
in the laboratory provided those SHN are not far away 
from /3-stability line. We have calculated the a decay 
half lives of SHN in quantum tunneling method with 
microscopic NN potential using Q-values from different 
mass formulae and compared them with the /3-decay 
and SF half lives to find the long lived SHN. 

The highlights of observations made in this work are 
summarised as follows: 

(i) Among all three mass formulae. Q^-values used 
from MMM model (Qm) 0, IM in the present 
method, reproduces the observed data reasonably well 
(see fig.l), but non-availability of Q^-values in more neu- 
tron rich side limits its usage. Therefore, for higher Z 
region, the mass formula of KUTY, which extends up to 
Z=130, has been considered. 

(ii) Although 266 Rf \q2 nucleus has relatively longer T a 
half life (~ 4.1 x 10 6 s —47.5 days using Q M ) but it is 
unstable against SF with Tsf SM- 23s only (Fig. 2). 

(iii) The mass formula of KUTY predicts a decay 
life time of the 290 Sgis4 nucleus to — 10 s yrs whereas 
Tsf SM~ 10 6 yrs makes the lifetime (~ 10 6 yrs) of this 
nucleus very long but still smaller than the age of the 
earth (~ 4.5 x 10 9 yrs) by three orders of magnitude. This 
nucleus is either /3-stable or might have very large Tjg ac- 
cording to the calculations of ref. [63| . 

(iv) The larger deviations between calculated and ex- 
perimentally measured a decay half lives are observed 
in case of only few nuclei such as 277 112 which may be 
due to higher minimum orbital angular momenta carried 
away by a particles for spin-parity conservation. Inaccu- 
racy in the measurement of T a of 277 112 nucleus due to 
very low count rate also may not be ruled out. 

(v) Usin g th e formulation based on liquid drop model 
in Ref. [641, |65| SF half lives calculation have also been 
done for higher Z elements in this work. Results shown 
in plots of figs. 4-5 do not match with SF half lives cal- 
culation in a microscopic approach of Ref.[6l], [62j for 
both neutron rich and neutron deficient isotopes of heav- 
iest elements with Z=104-120. The half life value using 
this phenomenological prescription also contradicts the 
observed a-decay from 283 112i7i nucleus with measured 
T a ~ 6.9s. 

(vi) It is evident from the present T a calculations that 
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the effect of the deformed neutron shell closure at N=162 
will be insignificant beyond Z=108 as T a goes on decreas- 
ing with increasing atomic number. For N=162 isotopes 
of elements having Z=110, 112, values of Tsf are 9.8m, 
0.63s respectively and T a using Qm are of the order of 
milliseconds (~ 1ms) and microseconds (~ 93/xs) (see 
plots of Fig. 3) respectively i.e. the values of T a go on 
decreasing more rapidly than the corresponding SF half 
lives with increasing atomic number. 

(vii) Calculated T a using Qkuty predicts almost 
/3-stable long lived SHN around 294 110i 84 , 296 H2i 84 , 
298 114i 84 with T Q of the order of - 311yrs, - 3.10yrs, ~ 
17 'days respectively which are much less than their Tsf 
(~ 4.48 x 10 4 yrs,~ 3.09 x 10 5 yrs, ~ 4.38 x I0 5 yrs re- 
spectively) values. Hence the dominant decay mode of 
the above nuclei and their immediate neighbours is ex- 
pected to have a decay mode. T a value of 293 110i83 is 
about 352 years which is slightly greater than that for 
294 110is4 nucleus. The SF half life of this nucleus is 
not found in Ref . [6ll. [62j]. For 292 108i 8 4 nucleus since 
T a (~ 9.6 x 10 yrs using Qkuty) value is comparable 
to its Tsf (~ 3.2 x 10 4 yrs), this nucleus is one of the 
possible members of stability island. The exact values 
of Tg of 293 110i83 and 292 108i84 nuclei are not shown in 
ref. [63] but predicted to be large. 

(viii) Using Qkuty values in the present calculation 
shows longer T a -values for neutron rich (N > 190) iso- 
topes of Z=116, 118 and 120 indicating a possible neu- 
tron shell closure next to N=184 might occur. On the 
contrary, calculated SF half lives (Fig. 3) of Ref. [HI, |62| 
show a trend of lowering of Tsf (< lms to lfis) for neu- 
tron rich isotopes of those elements which indicates the 
higher probability of SF of such SHN in this region. How- 
ever, more accurate determination of fission barrier and 
their corresponding half lives are essential to predict long 
lived SHN in the region of very high atomic number. 



(ix) It may be pointed out that the calculation of T a 
is very sensitive to Q^-values, and none of the mass 
formula used here cover the entire mass range with 
extreme accuracy. Therefore, a better mass estimate 
covering the wide range of superheavy masses with a 
good accuracy is necessary. 

In summary, we find that the possibility of existence 
of SHN above Z— 114 with considerable life time is very 
low. Although Z=120, 124, 126 with N=184 might form 
spherical-doubly-magic nuclei and survive fission [79l |. 
they would undergo a decay within microseconds. A 
small "island/peninsula" might survive fission and f3- 
decay but undergo a decay in the region Z=106-108, N 
~ 160-164. Interestingly, in this region the /3-stable SHN 
Z=106, N=162 has the highest a decay half life ~ 3.2hrs 
(Fig. 2, using Qm) that is much greater than the recently 
discovered dcformcd-doubly-magic SHN 270 Hs (measured 
T a ~ 22 sees). Thus a search for this long-lived SHN 
268 Sgi62 can be pursued. Similarly, the nucleus with 
Z=110, N=183 appears to be near the center of a pos- 
sible "magic island" (Z=104-116, N - 176-186) with a 
decay half life ~ 352yrs (Fig. 3, using Qkuty) which 
is greater than that of the doubly-magic SHN Z=114, 
N=184 (T a ~ 17days). Since the SHN 2m Sg 184 has 
T a and Tsf values ~ 10 s yrs (Fig. 2, using Qkuty), 
and ~ 10 6 yrs respectively, it might have longer life time 
in comparison to other superheavies. However, for both 
293 .Dsi83 and 290 S<7i84 nuclei, /3-decay might be another 
possible decay mode with large Tp values. Only future 
experiments can confirm this. Finally, the experimental 
investigations to detect the a-cascade can be pursued on 
294 110i84, 293 H0i83, 296 H2i84 and 298 114 184 nuclei which 
are expected to decay predominantly through a particle 
emission. 
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FIG. 1: "(Color online)" Plots of a decay half life [logio(Ti/2 / sec)] in logarithmic scale versus proton number Z for different 
mass number A (indicated on top of each coloumn) using zero angular momentum transfer (i = 0). (a) bar coded columns are 
theoretical half lives (T a M3Y Q-MS) in WKB frame work with DDM3Y interaction and [Q^ s ] from Myers-Swiatecki mass 
formula, (b) columns filled with dots (T a M3Y Q-K) are in the same framework but with [Q^, f/Ty ] from Koura-Tachibana- 
Uno-Yamada mass estimates, (c) solid columns are theoretical half lives (T a M3Y Q-M) in WKB frame work with DDM3Y 
interaction and [Qfl] from Muntian-Patyk-Hofmann-Sobiczewski mass formula, (d) hollow columns are experimental a decay 
half lives (T a Expt). Experimental errors are given in Table I. 
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FIG. 2: "(Color online)" Variation of a decay and fission half-lives with neutron number for elements (a) Z=102, (b) Z=104, (c) 
Z=106, (d) Z=108 are shown. For all plots the following symbols are used: Dash-dotted line (T a M3Y Q-K) and continuous line 
with square symbols (T a M3Y Q-M) represent a decay half-lives calculation using Q-values from KUTY (Q-K) and Muntian 
et al. (Q-M) respectively in this work. Triangle symbol (T a SM) represents a-decay half-lives predicted within a microscopic 
framework [oil. |62|]. Hexagon symbol (T a Expt) represents measured a decay half- lives. Solid circle (Tsf SM) represents fission 
half-lives predicted by microscopic calculation. Diamond symbol (Tsf expt) represents measured fission half-lives for some 
nuclei. Line-inverted triangle (Tp) shows /3 decay half-lives predicted in ref. [631 ]. 
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FIG. 3: "(Color online)" Same as Fig.2 for elements (a) Z=110, (b) Z=112, (c) Z=114, (d) Z=116, (e) Z=118, (f) Z=120 
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FIG. 4: "(Color online)" Plots for variation of phenomenological [64], [65j] and microscopically [oil . |62| calculated fission half 
lives with neutron numbers for (a) Z=102, (b) Z=104, (c) Z=106, (d) Z=108, (e) Z=110, (f) Z=112. The corresponding a 
decay half-lives from the present calculation are also shown for comparison. Continuous line with solid circle (T a M3Y Q-M) 
represents a decay half-lives predicted by this work using Q value from Muntian et al. 7, 68, 69]. Continuous line with solid 
triangle (Tsf Ph) represents spontaneous fission half-lives predicted by phenomenological calculation. Solid dark circle (Tsf SM) 
shows the spontaneous fission half-lives predicted by microscopic calculation. 




FIG. 5: "(Color online)" Same as fig.4 for (a) Z=114, (b) Z=116, (c) Z=118, (d) Z=120 
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FIG. 6: "(Color online)" Variation of a decay half-lives (T a ) predicted by present calculation with neutron number for (a) 
odd-Z and (b) even-Z elements from Z=102-120. Each graphs show the abrupt reduction of T a around N=162 as the possible 
signature of shell effect. Similar reductions of T a around N=184 are also shown by the elements for which Q a values are 
available from MMM calculation. 



